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Abstract 

Using a recent reformulation of the analysis of nuclear parity- violation 
(PV) using the framework of effective field theory (EFT), we show how 
predictions for parity-violating observables in low energy light hadronic 
systems can be understood in an analytic fashion. It is hoped that such 
an analytic approach may encourage additional experimental work as well 
as add to the understanding of such parity-violating phenomena, which 
is all too often obscured by its description in terms of numerical results 
obtained from complex two-body potential codes. 
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1 Introduction 



I never had the pleasure of meeting Dubravko Tadic , which is a shame, since 
he and I worked on many paraUel subjects over the years. An example of this is 
his recent work on hypernuclear decaypQ as well as his early papers on what is 
usually called nuclear parity violation j^j. It is the latter which I wish to focus 
on in this paper, which I dedicate to Dubravko's memory. 

The cornerstone of traditional nuclear physics is the study of nuclear forces 
and, over the years, phenomenological forms of the nuclear potential have be- 
come increasingly sophisticated. In the nucleon-nucleon (NN) system, where 
data abound, the present state of the art is indicated, for example, by phe- 
nomenological potentials such as AV18 that are able to fit phase shifts in the 
energy region from threshold to 350 MeV in terms of ^ 40 parameters. Progress 
has also been made in the description of few-nucleon systems j3| . At the same 
time, in recent years a new technique — effective field theory (EFT) — has been 
used in order to attack this problem using the symmetries of QCD jlj. In 
this approach the nuclear interaction is separated into long- and short-distance 
components. In its original formulation designed for processes with typical 
momenta comparable to the pion mass, Q ~ m^, the long-distance component is 
described fully quantum mechanically in terms of pion exchange, while the short- 
distance piece is described in terms of a small number of phenomenologically- 
determined contact couplings. The resulting potential [HUZ] is approaching JHUH] 
the degree of accuracy of purely-phenomenological potentials. Even higher pre- 
cision can be achieved at lower momenta, where all interations can be taken as 
short-ranged, as has been demonstrated not only in the NN system |10l 111) , 
but also in the three- nucleon system Precise — ~ 1% — values have 

been generated also for low-energy, astrophysically-important cross sections for 
reactions such asn + p^d + j |14| and p + p d-l-e"^-!- ^c|15|. However, 
besides providing reliable values for such quantities, the use of EFT techniques 
allows for the a realistic estimation of the size of possible corrections. 

Over the past nearly half century there has also developed a series of mea- 
surements attempting to illuminate the parity- violating (PV) nuclear interac- 
tion. Indeed the first experimental paper of which I am aware was that of Tanner 
in 1957 jl6) . shortly after the experimental confirmation of parity violation in 
nuclear beta decay by Wu et al. [El- Following seminal theoretical work by 
Michel in 1964 [TH| and that of other authors in the late 1960's [HOniEI], the 
results of such experiments have generally been analyzed in terms of a meson- 
exchange picture, and in 1980 the work of Desplanques, Donoghue, and Holstein 
(DDH) developed a comprehensive and general meson-exchange framework for 
the analysis of such interactions in terms of seven parameters representing weak 
parity-violating meson- nucleon couplings |22|. The DDH interaction has be- 
come the standard setting by which hadronic and nuclear PV processes are now 
analyzed theoretically. 

It is important to observe, however, that the DDH framework is, at heart, 
a model based on a meson-exchange picture. Provided one is interested primar- 
ily in near-threshold phenomena, use of a model is unnecessary, and one can 
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instead represent the PV nuclear interaction in a model-independent effective- 
field-theoretic fashion, as recently developed by Zhu et al.|22I- In this approach, 
the low energy PV NN interaction is entirely short-ranged, and the most general 
potential depends at leading order on 11 independent operators parameterized 
by a set of 11 a priori unknown low-energy constants (LEC's). When applied to 
low-energy (i^cm < 50 MeV) two-nucleon PV observables, however, such as the 
neutron spin asymmetry in the capture reaction n+p — s- d-l-7, the 11 operators 
reduce to a set of five independent PV amplitudes which may be determined by 
an appropriate set of measurements, as described in and an experimental 
program which should result in the determination of these couplings is underway. 
This is an important goal, since such interactions are interesting not only in their 
own right but also as background effects entering atomic PV measurements |22 
as well as experiments that use parity violation in electromagnetic interactions 
in order to probe nuclear structure |25j. 

Completion of such a low-energy program would serve at least three addi- 
tional purposes: 

i) First, it would provide particle theorists with a set of five benchmark num- 
bers which are in principle explainable from first principles. This situation 
would be analogous to what one encounters in chiral perturbation theory 
for pseudoscalars, where the experimental determination of the ten LEC's 
appearing in the 0{p'^) Lagrangian presents a challenge to hadron struc- 
ture theory. While many of the 0{p'^) LEC's arc saturated by t-channel 
exchange of vector mesons, it is not clear a priori that the analogous PV 
NN constants are similarly saturated (as assumed implicitly in the DDH 
model) . 

ii) Moreover, analysis of the PV NN LEC's involves the interplay of weak 
and strong interactions in the strangeness conserving sector. A similar 
situation occurs in AS = 1 hadronic weak interactions, and the interplay 
of strong and weak interactions in this case are both subtle and only 
partially understood, as evidenced, e.g., by the well-known the AI = 1/2 
rule enigma. The additional information in the AS = sector provided 
by a well-defined set of experimental numbers would undoubtedly shed 
light on this fundamental problem. 

iii) Finally, the information derived from the low-energy nuclear PV program 
would also provide a starting point for a reanalysis of PV effects in many- 
body systems. Until now, one has attempted to use PV observables ob- 
tained from both few- and many-body systems in order to determine the 
seven PV meson-nucleon couplings entering the DDH potential, and sev- 
eral inconsistencies have emerged. The most blatant is the vastly different 
value for obtained from the PV 7-decays of ^*F, ^^F and from the com- 
bination of the pp asymmetry and the cesium anapole moment j24|. The 
origin of this clash could be due to any one of a number of factors. Using 
the operator constraints derived from the few-body program as input into 
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the nuclear analysis could help clarify the situation. It may be, for exam- 
ple, that the remaining combinations of operators not constrained by the 
few-body program play a more significant role in nuclei than implicitly 
assumed by the DDH framework. Alternatively, truncation of the model 
space in shell model treatments of the cesium anapole moment may be 
the culprit. In any case, approaching the nuclear problem from a more 
systematic perspective and drawing upon the results of few-body studies 
would undoubtedly represent an advance for the field. 

The purpose of the present paper is not, however, to make the case for the 
effective field theory program — this has already been undertaken in [23]. Also, 
it is not our purpose to review the subject of hadronic parity violation — indeed 
there exist a number of comprehensive recent reviews of this subject[2SlEZIEll- 
However, although the basic ideas of the physics are clearly set out in these 
works, because the NN interaction is generally represented in terms of a some- 
what forbidding two-body interaction, any calculations which are done involve 
state of the art potentials and are somewhat mysterious except to those priests 
who preach this art. Rather, in this paper, we wish to argue that this need not 
be the case. Below we eschew a high precision but complex nuclear wavefunc- 
tion approach in favor of a simple analytic treatment which captures the fiavor 
of the subject without the complications associated with a more rigorous cal- 
culation. We show that, provided that one is working in the low energy region, 
one can use a simple effective interaction approach to the PV NN interaction 
wherein the the parity violating NN interaction is described in terms of just five 
real numbers, which characterize S-P wave mixing in the spin singlet and triplet 
channels, and the experimental and theoretical implications can be extracted 
within a basic effective interaction technique, wherein the nucleon interactions 
are represented by short range potentials. This is justified at low energy because 
the scales indicated by the scattering lengths — ~ — 20 fm, at ^ 5 fm — are 
both much larger than the ~ 1 fm range of the nucleon-nucleon strong inter- 
action. Of course, precision analysis should still be done with the best and 
most powerful contemporary wavefunctions such as the Argonne V18 or Bonn 
potentials. Nevertheless, for a simple introduction to the field, we feel that the 
elementary discussion given below is didactically and motivationally useful. In 
the next section then we present a brief review of the standard DDH formalism, 
since this is the basis of most analysis, as well as the EFT picture in which we 
shall work. Then in the following section we show how the basic physics of the 
NN system can be elicited in a simple analytic fashion, focusing in particular 
on the deuteron. With this as a basis we proceed to the parity-violating NN in- 
teraction and develop a simple analytic description of low energy PV processes. 
We summarize out findings in a brief concluding section 

2 Hadronic Parity Violation: Old and New 

The essential idea behind the conventional DDH framework relies on the fairly 
successful representation of the parity-conserving NN interaction in terms of a 
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Figure 1: Parity- violating NN potential generated by meson exchange. 



single meson-exchange approach. Of course, this technique requires the use of 
strong interaction couplings of the lightest vector and pseudoscalar mesons 



whose values are reasonably well determined. The DDH approach to the parity- 
violating weak interaction utilizes a similar meson-exchange picture, but now 
with one strong and one weak vertex — c/. Fig. 1. 

We require then an effective parity-violating NNM Hamiltonian in analogy 
to Eq. Q). The process is simplified somewhat by Barton's theorem, which 
requires that, in the CP-conserving limit, which we employ, exchange of neutral 
pseudoscalars is forbidden From general arguments, the effective Hamilto- 
nian for such interactions must take the form 



term is usually omitted, leaving parity-violating observables described in terms 
of just six constants. DDH attempted to evaluate such PV couplings using 
basic quark-model and symmetry techniques, but they encountered significant 
theoretical uncertainties. For this reason their results were presented in terms 
of an allowable range for each, accompanied by a "best value" representing their 
best guess for each coupling. These ranges and best values are listed in Table 
121 together with predictions generated by subsequent groups EH . 

Before making contact with experimental results, however, it is necessary to 
convert the NNM couplings generated above into an effective parity-violating 
NN potential. Inserting the strong and weak couplings, defined above into the 
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Table 1: Weak NNM couplings as calculated in Refs. |221|SI1E2- AH numbers 
are quoted in units of the "sum rule" value — 3.8 • 10~^. 



meson-exchange diagrams shown in Fig. 1 and taking the Fourier transform, 
one finds the DDH effective parity-violating NN potential 
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where Wi{r) = exp(— mir)/47rr is the usual Yukawa form, r = ja'i — X2\ is the 
separation between the two nucleons, and pi ~ — iV;. 

Nearly all experimental results involving nuclear parity violation have been 
analyzed using V^^^ for the past twenty-some years. At present, however, 
there appear to exist discrepancies between the values extracted for the various 
DDH couplings from experiment. In particular, the values of and /i^ ex- 
tracted from pp scattering and the 7 decay of ^^F do not appear to agree with 
the corresponding values implied by the anapole moment of ^^'^Cs measured in 
atomic parity violation j33) . 

These inconsistencies suggest that the DDH framework may not, after all, 
adequately characterize the PV NN interaction and provides motivation for our 
reformulation using EFT. In this approach, the effective PV potential is entirely 
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short-ranged and has the co-ordinate space form 
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and fm{T^ is a function which 

i) is strongly peaked, with width ^ 1/m about r = 0, and 

ii) approaches (5^'^-'(f*) in the zero width — to — > 00 — hmit. 

A convenient form, for example, is the Yukawa-like function 

^2 



exp{—'mr) 



(6) 



where to is a mass chosen to reproduce the appropriate short range effects. 
Actually, for the purpose of carrying out actual calculations, one could just as 
easily use the momentum-space form of Vg^ , thereby avoiding the use of fm{'r) 
altogether. Nevertheless, the form of Eq. 0]is useful when comparing with the 
DDH potential. For example, we observe that the same set of spin-space and 
isospin structures appear in both Vf^j and the vector-meson exchange terms in 
^DDHi though the relationship between the various coefficients in V^^^ is more 
general. In particular, the DDH model is tantamount to assuming 
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(7) 
(8) 



C3 _ C4 _ C5 
C3 C4 C5 

and taking to ~ mp,mui, assumptions which may not be physically realistic. 
Nevertheless, if this ansatz is made, the EFT and DDH results coincide provided 
the identifications 
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are made '23'. 

Before beginning our analysis of PV NN scattering, however, it is important 
to review the analogous PC NN scattering case, since it is more familiar and it 
is a useful arena wherein to compare conventional and effective field theoretic 
methods. 



3 Parity Conserving NN Scattering 

We begin our discussion with a brief review of conventional scattering theory [531 ■ 
In the usual partial wave expansion, we can write the scattering amplitude as 

f{e) = ^(2£+ l)a,(fc)P,(cos0) (9) 
e 

where a^(fc) has the form 

°''^''r'""-'°^'^'- tco.^-a^ (10) 

3.1 Conventional Analysis 

Working in the usual potential model approach, a general expression for the 
scattering phase shift Si{k) is 34 



smSi{k)^-k dr'r'jeikr')2mrV{r')ue,k{r') (11) 
Jo 

where mr is the reduced mass and 

ut^k{r) = r cos Si{k)ji{kr) + kr dr'r' ji{kr')ni{kr)ui^k{r')2mrV{r') 



+ kr dr'r'je{kr)n{{kr')ue^k{r')'2mrV{r') (12) 



7 



is the scattering wavefunction. At low energies one can characterize the analytic 
function cot S{k) via an effective range expansion |35| 

k^''+^cotSe{k) = -^ + ^r,k^ + ... (13) 
Then from Eq. Illl we can identify the scattering length as 

= r.o. ,^,n9 / dr' {r' 2mrV {r' ) + ©(F^) (14) 



[(2^+l)!!P7o 

For simplicity, we consider only S-wave interactions. Then for neutron-proton 
interactions, for example, one finds 

ag = -23.715 ± 0.015 fm, = 2.73 ±0.03 fm 

al = 5.423 ± 0.005 fm, 1.73 ± 0.02 fm (15) 

for scattering in the spin-singlet and spin-triplet channels respectively. The 
existence of a bound state Eb = — 7^/2mr is indicated by the presence of a pole 
along the positive imaginary fc-axis — i.e. 7 > under the analytic continuation 
k ij — 

- + -7 = (16) 

We see from Eq. El that there is no bound state in the np spin-singlet channel, 
but in the spin-triplet system there exists a solution 



1- /l-^ 

^ ^ V <_ ^ 45 7 ]^gY^ ■ ^ _2.23 MeV (17) 

corresponding to the deuteron. 

As a specific example, suppose we utilize a simple square well potential to 
describe the interaction 

For S-wave scattering the wavefunction in the interior and exterior regions can 
then be written as 

^^o(ifr) r<R 
i N' {jo{kr) cos So ~ no{kr) sin Sq) r>R ^ ^ 

where jo, no are spherical harmonics and the interior, exterior wavenumbers 
are given by A: = ^/2mrE, K — ^2mr{E + Vo) respectively. The connection 
between the two forms can be made by matching logarithmic derivatives at the 
boundary, which yields 



Making the effective range expansion — Eq^] — we find an expression for the 
scattering length 
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where Kq ~ \/2mrVo 



(21) 



Note that for weak potentials — KqR << 1 — this form agrees with the general 
result Eq. [HI— 



ao 



dr'r' 2mrV{r') = 



2m 



r d3 



R'Vo+OiVo') 



(22) 



3.2 Coulomb Effects 

When Coulomb interactions are included the analysis becomes somewhat more 
challenging. Suppose first that only same charge {e.g., proton-proton) scattering 
is considered and that, for simplicity, we describe the interaction in terms of a 
potential of the form 

U[r) r <R 
^ r> R 



Vir) 



(23) 



i.e. a strong attraction — U{r) — at short distances, in order to mimic the strong 
interaction, and the repulsive Coulomb potential — a/r — at large distance, where 
a ~ 1/137 is the fine structure constant. The analysis of the scattering then 
proceeds as above but with the replacement of the exterior spherical Bessel 
functions by corresponding Coulomb wavefunctions F^^ , Gq 



(24) 



whose explicit form can be found in reference |3()| . For our purposes we require 
only the form of these functions in the limit kr « 1 — 
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Here 7^ = 0.577215.. is the Euler constant, 

2'KX 



exp(27ra;) — 1 



(25) 



(26) 



is the usual Coulombic enhancement factor, as — l/rrira is the Bohr radius, 
?7+(fc) = l/2fcas, and 



h{v+{k)) = ReH{irj+{k)) = rjl{k) ^ 



n(n2 +7^2 (fc)) 



lnv+{k)-jE (27) 
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where H(x) is the analytic function 



H{x) - iPix) + ^ - Inix) (28) 



Equating interior and exterior logarithmic derivatives we find 
KF{KR) 



cosSoF+'{R) - sin(5oG(['(i?) 



cos 5aF+{R) - sin SaG^{R) 
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(29) 

Since R « as Eq. [23 can be written in the form 
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The scattering length ac in the presence of the Coulomb interaction is conven- 
tionally defined as|37| 

kcotSoC^iv+ik)) + —h{Tj+{k)) = -— + ... (31) 
as ac 

so that we have the relation 

-- = ---- (ln^ + l-27s) (32) 
ao ac CLB R 

between the experimental scattering length — ac — and that which would exist 
in the absence of the Coulomb interaction — ag. 

As an aside we note that, strictly speaking, ap is not itself an observable 
since the Coulomb interaction cannot be turned off. However, in the case of 
the pp interaction isospin invariance requires a^ — a™ so that one has the 
prediction 

While this is a model dependent result, Jackson and Blatt have shown, by 
treating the interior Coulomb interaction perturbatively, that a version of this 
result with 1 — 2^e 0.824 — 2^e is approximately valid for a wide range 
of strong interaction potentials and the correction indicated in Eq. is 
essential in restoring agreement between the widely discrepant — a™ = —18.8 
fm vs. aj^ — —7.82 fm — values obtained experimentally. 

Returning to the problem at hand, the experimental scattering amplitude 
can then be written as 



where uo = argr(l — irj+{k)) is the Coulomb phase 



2*'^«C2(?7+(fc)) 
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3.3 Effective Field Theory Analysis 

Identical results may be obtained using effective field theory (EFT) methods 
and in many ways the derivation is clearer and more intuitive |38|. The basic 
idea here is that since we are only interested in interactions at very low energy, 
a scattering length description is quite adequate. From Eq. 122 we see that, at 
least for weak potentials, the scattering length has a natural representation in 
terms of the momentum space potential V{p — 0) — 



(35) 



and it is thus natural to perform our analysis using a simple contact interation. 
First consider the situation that we have two particles A,B interacting only via 
a local strong interaction, so that the effective Lagrangian can be written as 



(36) 



The T-matrix is then given in terms of the multiple scattering series shown in 
Figure 1 



Tf,{k) = -— /(fc) = Co + C^Goik) + C^Glik) + ... = - ^° 
rrir 1 



CoGo(fc) 



(37) 



where G'o(fc) is the amplitude for particles A, B to travel from zero separation 
to zero separation — i.e the propagator Dp{k;r' = 0, r = 0) — 



G'o(fc) = lim / -77r^-r2 -2 = / 7:7^ 



2m„ 



Equivalently Tfi{k) satisfies a Lippman-Schwinger equation 

T/,(fc) = Co + CoGo(fc)ry,(fc). 
whose solution is given in Eq. 1371 



(38) 



(39) 




Figure 2: The multiple scattering series. 



+ 



The complication here is that the function Go{k) is divergent and must be 
defined via some sort of regularization. There are a number of ways by which to 
do this, but perhaps the simplest is to use cL cutoff regulcLrization with kmax — 
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which simply ehminates the high momentum components of the wavcfunction 
completely. Then 

Go{k) = -^{^-^+^k) (40) 

(Other regularization schemes are similar. For example, one could subtract at 
an unphysical momentum point, as proposed by Gegelia[5^ 

^ /,\ f , 2mr 2to,. , nir , .,, 

which has been shown by Mehen and Stewart (lU] to be equivalent to the power 
divergence subtraction (PDS) scheme proposed by Kaplan, Savage and Wise. [55]') 
In any case, the would-be linear divergence is, of course, cancelled by introduc- 
tion of a counterterm accounting for the omitted high energy component of the 
theory, which renormahzes Cq to Co(/x). That Co(^) should be a function of the 
cutoff is clear because by varying the cutoff energy we are varying the amount 
of higher energy physics which we are including in our effective description. The 
scattering amplitude then becomes 

Comparing with Eq. ^| wc identify the scattering length as 

ao mrCo{ii) TT 

Of course, since oq is a physical observable, it is cutoff independent, so that the 
fj, dependence of l/Co(/i) is cancelled by the cutoff dependence in the Green's 
function. 



3.4 Coulomb Effects in EFT 

More interesting is the case where we restore the Coulomb interaction between 
the particles. The derivatives in Eq. 1361 then become covariant and the bubble 
sum is evaluated with static photon exchanges between each of the lines — each 
bubble is replaced by one involving a sum of zero, one, two, etc. Coulomb 
interactions, as shown in Figure 2. 





+ 




+ 




+ 



Figure 3: The Coulomb corrected bubble. 



The net result in the case of same charge scattering is the replacement of 
the free propagator by its Coulomb analog 



Go(fc) -> G+(fc) 



lim 

r', r— >0 



(27r)3 _ + ^e 
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d^s 2mrCHv+{s)) ^^^^ 



(27r)3 k^-s^+ie 



where 

'iP+{r)^C{r]+{s))e"'°e'^-^iFi{-iri+{s),l,isr~is-r) (45) 

is the outgoing Coulomb wavefunction for repulsive Coulomb scattering. (41) Also 
in the initial and final states the influence of static photon exchanges must be 
included to all orders, which produces the factor C^(27r77+(fc)) exp(2icro). Thus 
the repulsive Coulomb scattering amplitude becomes 

f+(k) = ™'' ^o'^^(^+(fe))^xP^^<^o C4g^ 
■^'^^ ' 2tt l-CoG+{k) ^ ' 

The momentum integration in Eq. 1441 can be performed as before using cutoff 
regularization, yielding 



r ( 2fl 1 



GUk)--^[^ + ^[H{^vAk))-\n^^-C } (47) 
where C = In 27r — 7. We have then 



C2(?7+(fc))e2^'^o 



- ^ Hv+ik) - In - C] - ^kC^iv+ik)) 
Comparing with Eq. |^ we identify the Coulomb scattering length as 



(48) 



-^ = -^ + ^(ln^ + C) (49) 
ac ao as tt 

which matches nicely with Eq. if a reasonable cutoff /i is 
employed. The scattering amplitude then has the simple form 

in agreement with Eq. 1^ 

Before moving to our ultimate goal, which is the parity violating sector, it is 
useful to spend some additional time focusing on the deuteron state, since this 
will be used in our forthcoming PV analysis and provides a useful calibration 
of the precision of our approach. 



4 The Deuteron 

Fermi was fond of asking the question "Where's the hydrogen atom for this 
problem?" meaning what is the simple model that elucidates the basic physics 
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of a given systemg3]? In the case of nuclear structure, the answer is clearly the 
deuteron, and it is essential to have a good understanding of this simplest of 
nuclear systems at both the qualitative and quantitative levels. The basic static 
properties which we shall try to understand are indicated in Table 2. Thus, 
for example, from the feature that the deuteron carries unit spin with positive 
parity, angular momentum arguments demand that it be constructed from a 
combination of S- and D-wave components (a P-wave piece is forbidden from 
parity considerations — more about that later). Thus the wavcfunction can be 
written in the form 

'^'^^^^ " 7^ (^"'i(^) + ■^Wd(r)Op„^ xt (51) 
where xt is the spin-triplet wavcfunction and 

is the tensor operator. Here Ud{r),Wd{r) represent the S-wave, D-wave compo- 
nents of the deuteron wavcfunction, respectively. We note that 

Opn\n> = (cos2 0-i)| tT>+sin2ee^'^Ui> 



Using 



?sin0e'^(| ti>+UT>) (52) 

dn. . 1 



Att ' ' 3 



S,, (53) 



we find the normalization condition 



< M^d >= y dr J dn [uiir) 
9 / 1 

-w}{r) {cos^e- -)2 + 2cos^6lsin2 6i + sin'*( 
8 \ 3 

/ 9 2 1 

dr U2(r) + -u;^(r)(l-- + -) 

dr{ulir)+w-'Ar)) (54) 

In lowest order we can neglect the D-wave component Wd{r). Then, in the region 
outside the range ro of the NN interaction we must have 

r>ro (-^^ + g)-«-0 (55) 

where 7 = 45.3 MeV is the deuteron binding momentum defined above. The 
solution to Eq. Ellis given by 

r > ro Ud{r) ^ e"^*" (56) 
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Binding Energy Eb 2.223 MeV 

Spin-parity 1+ 

Isospin T 

Magnetic Dipole Moment /id 0.856/iAr 

Electric Quadrupole Moment Qd 0.286 efm^ 

Charge Radius ^Jr^ ^ 2 fm 

Table 2: Static properties of the deuteron. 



However, since I/7 ~ 4.3 fm >> ro ^ 1 fm, it is a reasonable lowest order 
assumption to assume the that the deuteron wavefunction has the form Eq. 1561 
everywhere, so that wc may take 

Of course, we also must consider scattering states. In this case the asymp- 
totic wavefunctions of the ^Si and ^Sq states must be of the form 

^(+) (r) '™ 4^ sin(/tr + 5{k)) = + —t{k) (58) 

kr kr r 

with 

t{k) = ie'*(^) sin6{k) 
k 

being the partial wave transition amplitude. At very low energy we may use 
the simple effective range approximation defined above 

A:ctn(5(fc) ~ -i (59) 



to write 



Then we have the representation 

^(+) (,) f:!^ sin(fcr + 5{k)) = ^ + — { (61) 
kr kr r \ — ik I 

and, comparing with a Green's function solution to the Schrodinger equation 

Ti.r /• ik\r-r'\ 
V,(+)(f) = e"^-'- — J d\'y—^U {7^)^(7^) (62) 

we see that the potential at low energy can be represented via the simple local 
potential 

4-Tr 

U{7^ ^ -ad^r) (63) 
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which is sometimes called the zero-range approximation (ZRA) and is equivalent 
to the contact potential used in the EFT approach. 

The relation between the scattering and bound state descriptions can be ob- 
tained by using the feature that the dcuteron wavefunction must be orthogonal 
to its '^Si counterpart. This condition reads 



= y d^rip^i (r ) Vt (r) = ^ '^^^'^'^ {I ^™ '^^ 



(fc) 



which requires that 7 = 1/ot. This necessity is also clear from the already 
mentioned feature that the deuteron represents a pole in tt{k) in the limit as 
k — > 17 — i.e., — 1/04+7 = 0- Since i ^ 4.3 fm, this equality holds to within 20% 
or so and indicates the precision of our approximation. In spite of this roughness, 
there is much which can be learned from this simple analytic approach. 
We begin with the charge radius, which is defined via 

< rl >=< 4>+l d\^r^\Mr)\^ (65) 

Note here that we have included the finite size of the proton, since it is compa- 
rable to the deuteron size and have scaled the wavefunction contribution by a 
factor of four since Tp = ^f. Performing the integration, we have 

J d^r^r^\Mr)\^ = w drr^u\r) = (66) 
and, since < >~ 0.65 fm^ we find 



<rj> = ^0.65 + ^ fm ~ 1.8 fm (67) 

which is about 10% too low and again indicates the roughness of our approxi- 
mation. 

Now consider the magnetic moment, for which the relevant operator is 

-* 6 C 

= i^v{Sp - Sn) + (MS - \){5p + CT„)^ (68) 

where J = L+i((Tp-|-(T„) is the total angular momentum, — A*p^Mn — 4.70, 
A*s = IJ'p+IJ'u = 0.88 are the isovector, isoscalar moments, and the "extra" factor 
of 1/2 associated with the orbital angular momentum comes from the obvious 
identity Lp = L/2. We find then 

= <V'd;i,i|M3|Vd;i,i>=^[i + (2Ms-i)<V'd;i,i|^3|V'd;i,i>] 
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e 
4M 

e 
2M 



Ms 



^(Ms - ^) / drwj{r) 



(69) 



In the lowest order approximation — neglecting the D-wave component of the 
deuteron — we find 

<l,l|M3|l,l>^Ms^ (70) 

and this prediction — fj,d = fj-s = — is in good agreement with the exper- 

imental value jU^^^ = 0.856/XAr- 

4.1 D-Wave Effects 

A second static observable is the quadrupole moment Qd, which is a measure 
of deuteron oblateness. In this case a pure S-wave picture predicts a spherical 
shape so that Qd = 0. Thus, in order to generate a quadrupole moment, we 
must introduce a D-wave piece of the wavefunction. Now, just as we related the 
£ = wavefunction to the np scattering in the spin triplet state, we can relate 
the D-wave component to the scattering amplitude provided we include spin. 
Thus, if we write the general scattering matrix consistent with time reversal and 
parity-conservation as|43j 



A4{k', k) = a + l3Sp ■ n3n ■ fi + p{3p + On) ■ n 

+ {k + X)(7p ■ fi^Sn • + [k — X)Sp ■ fi_(7„ • fl^ 



(71) 



where 



k±k' 



kxk' 



\k±k'\ \kxk'\ 
we can represent the asymptotic scattering wavefunction via 

^P{r)"^ e'^-'' + Mi-iV,k)^ (72) 
A useful alternative form for M can be found via the identity 

(jp ■ nan ■ n = (jp ■ (In — fJp ■ n+cjn ■ n+ — fJp ■ fi-ffn ■ fi- (73) 
and, using the deuteron spin vector S — \{(Jp + dn), it is easy to see that 

M{k', k) = -ot + -i? \c'S -kxk' + gi{S-{k + k')f + 92(8 ■ {k - k'))^] (74) 



where 



2pAP 
fc^ sin 6 ' 



51 



2fc2 C0S2 ^ 



52 = 



(k - /3 - X)AP 
2fc2 sin^ 
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Then, since in the ZRA 

k' -^iVS^if). 
we find the effective local potential 



(75) 



47r 
1 



(76) 



where 



4 

Sij — SiSj + SjSi — g'^ij 



(77) 



Using the Green's function representation — Eq. [721 the asymptotic form of the 
triplet scattering wavefunction becomes 



, 31+52 C Y7 Y7 



ikr 



-Xt 



(78) 



and, by continuing to the value k — > 17, we can represent the deuteron wave- 
function as 

27r 



Mr) 



(79) 



A little work shows that this can be written in the equivalent form 



2_ 

2ti 



1 I 91+92 



(P 1 d 
dr"^ r dr 



-e-^'Xt 
r 



1 I 91+92 



37 



-7 



(80) 



Here the asymptotic ratio of S- and D-wave amplitudes is an observable and is 
denoted by 

Ad _ V2(5i+52) 



A, 



(81) 



This quantity has been determined experimentally from elastic dp scattering 
and from neutron stripping reactions to be|44| 

r? = 0.0271 ± 0.0004, i.e., gi + 32 = 105 fm 

Defining the quadrupole operator^ 



^Notc that the factor of j arises from the identity Vp = 



18 



and using 

—firjfkfe = j^{5ijSke + SuSjk + SikSje) (82) 
47r 15 



/ 



we note 



37 , r 



3e 1 51+92)7 7 t , - - ^ 

^ T ' 15 — 2M2 — 277 I o-pi<T„j + apjani - -^dijap -anjXt 

X 4tt fire"2"''''(3 + 37r + 7^^) 
Jq 

1 e(gi + 32)7' t ^ ^ 2 - \ / 3 37 ^ 272 

= 5 2M5 ^* V"'^'''"'' " ■ ^* (2^ (2^ 

e(gi + g2)7 t / _L 2^ - - ^ 

= -^j^ Xt I CTpicmj + cTpjani - -^Oijap ■ an j Xt 



Then the quadrupole moment is found to be 

Qf =< 1, l\Qzz\^d; 1, 1 >= ^^^^ - 0.28 e W, (84) 
in good agreement with the experimental value 

Q^'^P = 0.286 

From its definition, we observe that the quadrupole moment would vanish for a 
spherical (purely S-wave) deuteron and that a positive value indicates a slight 
elongation along the spin axis. 

Note that in interpreting the meaning of the D-wave piece of the deuteron 
wavefunction one often sees things described in terms of the D-state probability 

/•oo 

Pd= drw'^ir) 
Jo 

/ + 7- + ^')' exp(-27r) (85) 

diverges while in reality the D-wave function Wd(r) must vanish as r ^ 0, it is 
clear that the connection between the asymptotic amplitude 77 and the D-state 
probability Pd must be model dependent. Nevertheless, the D-state piece is a 
small but important component of the deuteron wavefunction. As one indication 
of this, let's return to the magnetic moment calculation and insert the D-wave 
contribution. We find then 

3 1 

Md = Ms - 2 il^s - -^)Pd (86) 



However, since 
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If we insert the experimental value fid — 0.857 we find Pd — 0.04 which can 
now be used in other venues. However, it should be kept in mind that this 
analysis is only approximate, since we have neglected relativisitic corrections, 
meson exchange currents, etc. 

Of course, static properties represent only one type of probe of deuteron 
structure. Another is provided by the use of electromagnetic interactions, 
for which a well-studied case is photodisintegration — 7c? np — or radiative 
capture — np — > dj — which are related by time reversal invariance. 



5 Parity Conserving Electromagnetic Interaction: 

np ^ 

An important low energy probe of deuteron structure can be found within the 
electromagnetic transtion np d'y. Here the np scattering states include both 
spin-singlet and -triplet components and we must include a bound state — the 
deuteron. For simplicity, we represent the latter by purely the dominant S-wave 
component, which has the form 

Mr) = \f^-e-'''-, or = (87) 

1' 27r r 7 -I- q 



Since we are considering an electromagnetic transition at very low energy we 
can be content to include only the lowest — El, Ml, and E2 — mutipoles, which 
are described by the Hamiltonian |42| 



Here Sj is the photon momentum, /iy, fig = fip^ fin are the isoscalar, isovector 
magnetic moments, and we have used Siegert's theorem to convert the conven- 
tional p- A/M interaction into the El form given above. ^ The ± in front of the 
El operator depends upon whether the np — *■ ^7 or ^d — > np reaction is under 
consideration. The electromagnetic transition amplitude then can be written in 
the form 

Amp = x\[^-t^s^-{GMiv{dp-dn) + GMis{Sp + Sn)) 

) X^ (89) 

The leading parity-conserving transition at near-threshold energy is then the 
isovector Ml amplitude — Gmiv — which connects the '^S'l deuteron to the ^5*0 
scattering state of the np system. From Eq. |SSlwe identify 

^^"^ = ^ / d'r^^^*ikr)Mr) (90) 



■^Note that the factor of two (eight) in the El (E2) component arises from the obvious 
identity irl^. 
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Using the asymptotic form 

-iS 

the radial integral becomes 



ipi g{kr) = — — (sin A:r cos(5s + cos A:r sincJs) (91) 



d'^ripigj {kr)^d{r) — — -j^ — J (sin fcr cos 5^ + cos A;r sin 5s)e 

-(fccosJg + 7sin^s) (92) 



Since by energy conservation 



k{k^ + 7^) 



^2 + 7^ 



we can use the lowest order effective range values for the scattering phase shift 
to this result in the form 

^ ^ eMy\/8^e-^*"""'fe°-(l-7a«) .gg. 

AM^^l + k-'aj ^ ' 

Note here that the phase of the amplitude is required by the Fermi- Watson 
theorem, which follows from unitarity. 

The Ml cross section is then found by squaring and mutiplying by phase 
space. In the case of radiative capture this is found to be 



Here 

Y^-TrPtTT^ = i^^^^Tri(3 + ai-a2)e;xs^-(ap-a„)e%xs^-(ap-a„) 
= ^IGMiype^ X -e^ X = 2|GMiyp (95) 

yielding 

aMi{np^dj) = ^-^^2\GMiv\ = ^^(1 + ^^,2) O^) 

Putting in numbers we find that for an incident thermal neutrons with relative 

velocity \vrei \ = 2200 m/scc. the predicted cross section is about 300 mb which 
is about 10% smaller than the experimental value 

o-e^p = 334 ± 0.1 mb 
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(a) (b) 




The discrepancy is due to our omission of two-body effects (meson exchange 
currents) as shown by Riska and Brown 45 . 

In a corresponding EFT description of this process, we must calculate the 
diagrams shown in Figure 1. There is a subtlety here which should be noted. 
Strictly speaking, as shown by Kaplan, Savage, and Wise^Bl the symbol (g) in 
these diagrams should be interpreted as creation or annihilation of the deuteron 
with wavefunction rcnormalization 



followed by propagation via 



1 



M M 



However, since in lowest order we have 

/dS(i;)\"^ 



we find for the product 



(98) 



= (99) 

M ^ M ' 

which is the deuteron wavefunction in momentum space. Thus in our discus- 
sion below we shall use this substitution rather that writing the wavefunction 
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normalization times propagator product. From Figm'e 4a then we find 



Gl 



MIV 



espfly 

4M 



{2n) 



3 



(100) 



Since {q) — (27r)^(5^(fc — q) we have 



Ga 



AM 72 + A;2 
On the other hand from Figures 4b+c we find 



'■^MIV 



AM 



Cos 



Since 



this becomes 



AM 

eSQfly 

AM 
esony 
AM 



l-Co.Go(fc) 



l-Co.Go(fc) 

Cqs 

l-CosGo(k) 

Cqs 



_fq_ 

(2^)3 

1 



Goif=0,q)Mq) 



M M ^ 



(101) 

(102) 
(103) 



/Wry 



(27r)3 (g2+^2)(|i„|i+,,) 

Wry 



ry2 _|_ ^2 

y/Wn M 



l-CosGo{k)J 72 + fc2 4^ 



(Go(fc)-Go(^7)) 
(~i/c — 7) 



(104) 



Adding the two contributions we have 

eSofiy y/Wrf 



G 



MIV 



AM 72 + fc2 



—7 — ik 
^-ik 



esofiy ^/Wrf /I — 7a, 



4M 72 + ^2 \ 1 + ika. 



(105) 



which agrees completely with Eq. l93l obtained via conventional coordinate space 
procedures. Of course, we still have a ~ 10% discrepancy with the experimental 
cross section, which is handled by inclusion of a four-nucleon Ml counterterm 
connecting ^Si and ^Sq states — 



(106) 



where here the Pi represent relvant projection operators|47|. 

As the energy increases above threshold, it is necessary to include the cor- 
responding P-conserving El multipole — 



Amp = Gfiie^ • kxlxt 



(107) 
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In this case the matrix element involves the np ^P-wave final state and, neglect- 
ing final state interactions in this channel, the matrix element is given by 

GEik=^ J d\i^f*{r)rMr) (108) 

The radial integral can be found via 

- ik ■ f d^re-'^-^fipd{r) = 4r I (fire-'^'^^-e-'^'' 

J dX\^^^ \ 2tt J r 

Equivalently using EFT methods we have 

J (2^0^"^^ (9)v,-^d(g) = (fe2 + ^2)2 (110) 

In either case 

_-iesoV8n2 ..... 

- (fc2+^2)2 (111) 

and the corresponding cross section is 



a^r(np-d7) = / ^(^^ " " (^^ " ^)') 



(fc2 + ^2)4 

87rafc27 



|tT,edM3(A;2+72) 



(112) 



It is important to note that we can easily find the corresponding photodis- 
integration cross sections by multiplying by the appropriate phase space. For 
unpolarized photons we have 

- 3^2^/(0i:2-^(^o4-£)|Amp|2 

■E/^IW (113) 



Mk 

247rso , 

A-y 



Using the results obtained above for radiative capture — 

|Amp|2 = 8|GMiy|' + 3(fc - - {k-sf)\GEi? (114) 

we find the photodisintegration cross sections to be 

27ra (l-as7)V^fc7 , , . STrajk^ 

CTMI [jd np) = /,.2 I -.2^n , 7.2.2^ ' [l/d np) = 



3M2 (fc2 + ^2) (1 + ^2^2)' 3(^2+^2)3 

(115) 
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Although the leading electromagnetic physics is controlled, as we have seen, 
by the isovector Ml and El amplitudes, there exist small but measurable 
isoscalar Ml and E2 transitions |3S],|13,|Sn]- In the former case the transition is 
between the S-wave (D-wave) deuteron ground state and into the ^Si{^D) scat- 
tering state. The amplitude Gmis is small because of the smallness of /Ug — ^ 
and because of the orthogonality restriction. In the case of G_e2, the result 
is suppressed by the requirement for transfer of two units of angular momen- 
tum, so that the transition must be between S- and D-wave components of the 
wavefunction. 

We first evaluate the isoscalar Ml amplitude, which from Eq. IHHlis given by 

(c/. Eq. inni) 

e 1 

e I [°° 1 

= 2m'^^^ ~ 2-* / ^^("d(r)ut(r) - -Wd{r)wt{r)) 

e 1 3 f°° 

^ ~2m'-^^ ^ 2''2 X ^^^^^^^'^'^'^^^ 

where the last form was fouind using the orthogonality condition. In order 
to estimate the latter, we follow Danilov in assuming that, since the radial 
integral is short-distance dominated, the D-wave deuteron and scattering pieces 
are related by a simple constant which, using orthogonality, must be given 

by 

Wt{r)~~atJ — Wd{r) (117) 
V 7 

The matrix element then becomes 

GMi5-^(M5-^)^Pz.a. (118) 
Likewise, using Eq. |H21we note that 



3 1 



drr'^(ut{r)wd(r) -|- Wt{r)ud{r)) [ap ■ e^iT„ ■ + i 



71 ^j^p -^7] 

(119) 



SO that the corresponding E2 couphng is found to be 



poo I 

"^'^i^-l *--n«.M».W + ..,M...M)^-^^yjf (120) 

In order to detect these small components we can use the circular polariza- 
tion induced in the final state photon by an initially polarized neturon, which 
is found to be 

Gmis — Ge2 



^ „2 ^'-''^ (121) 
Gmiv 



25 



Putting in numbers we find 

" ^,v{l-la,)V' 15 atA-P ) 

~ -1.17x10^^-0.24x10"^ 

= -1.41 X IQ-^ (122) 
which is in reasonable agreement with the experimental value [K^ 

pexp ^ 5 ^ Q 3) ^ (^23) 

Having familiarized ourselves with the analytic techniques which are needed, 
we now move to our main subject, which is hadronic parity violation in the NN 
system. 



6 Parity-Violating NN Scattering 

For simplicity we begin again with a system of two nuclcons. Then the NN 
scattering-matrix can be written at low energies in the phenomenological form|51| 



Mik', k) = mt{k)Pi + ms{k)PQ 



(124) 



where 



Pi - J (3 + (71 -^2), Po = ^(l-'?l-f?2) 



are spin-triplet, -singlet spin projection operators and 

-at -a.s 



mtik) = — 



ikat ' 



-|- ikug 



(125) 



are the S-wave partial wave amplitudes in the lowest order effective range ap- 
proximation, keeping only the scattering lengths at, as- Here the scattering cross 
section is found via 

^ = TiM^M (126) 
ail 

so that at the lowest energy we have the familiar form 



dn 



1 + 



(127) 



The corresponding scattering wavefunctions are then given by 



4+)(r) 



Ak-r 



l\'f r ^ik\r—r'\ ^ ^ 



Ak-r 



Mi-iV,k)- 



Akr " 



(128) 



26 



where x is the spin function. In Born approximation we can represent the 
wavefunction in terms of an effective deha function potential 

UtA^) = + asPo)6^{r) (129) 

as can be confirmed by substitution into Eq. 11281 

6.1 Including the PV Interaction 

Following Danilov.|51|. we can introduce parity mixing into this simple represen- 
tation by generalizing the scattering amplitude to include P-violating structures. 
Up to laboratory energies of 50 MeV or so, we can omit all but S- and P-wave 
mixing, in which case there exist only five independent such amplitudes: 

i) dt(k) representing '^Si ^Pi mixing; 

ii) d°'^'^(A;) representing ^5*0 -^Pq niixing with A/ = 0, 1, 2 respectively; 

iii) Ct(fc) representing ^Si '^Pi mixing. 

After a little thought, it becomes clear then that the low energy scattering- 
matrix in the presence of parity violation can be written as 

M{k',k) = m,{k)Po + Ct{k){ai+a2)-{k' + k)^{Ti-T2), 

+ (ffi S2) ■ (k' + k) (P,dl{k) + i(ri + T,)^dl{k) + ^-Il^I^^^S^dlik)^ 
mt{k) + dt{k){ai ~ a^) ■ {k' + fc)] Pi (130) 

Note that since under spatial inversion — a — > (7,k,k' —k,~k' — each of the 
new pieces is P-odd, and since under time reversal — a —^ —(j,k,k' — > —k',—k 
the terms are each T-even. At very low energies the coefficients in the T-matrix 
become real and we define [ST] 

limms^tik) = as,t, lim ct{k), ds{k), dt{k) ^ ptttt, Xla^, XtUt (131) 

(The reason for factoring out the S-wave scattering length will be described 
presently.) The five real numbers pt, A^, At then completely characterize the low 
energy parity-violating interaction and can in principle be determined experi- 
mentally, as we shall discuss below. Alternatively, we can write things in terms 
of the equivalent notation 

Af = A^ + Ai + -i=A^ 



•^Note that there exists no singlet analog to the spin-triplet constant ct since the combina- 
tion (Ti -I- (T2 is proportional to the total spin operator and vanishes when operating on a spin 
singlet state. 
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Ar = A°-Ai + -^A^ (132) 

We can also represent this interaction in terms of a simple effective NN 
potential. Integrating by parts, we have 

/pik\r — r'\ _ pikr 

= + (133) 

which represents the parity-violating admixture to the the scattering wavcfunc- 
tion in terms of an S-wave admixture to the scattering P-wave state — ~ a ■ k^-y- 

plus a P-wave admixture the scattering S-state — <^ —ia ■ V^^. We see then 
that the scattering wave function can be described via 

C/(r) = ^ [{acS-'{r) + XtUti^i - aa) • {-»V, S'^ir)}) P, 

+ as6^f)Po + Ptat{ai + az) • {-«V, <53(f)}i(ri - ra)^ 

+ (ai - a2) ■ {-iW, SHf^}a, (PoA° + ^(n + r,)X + ^^^^l^^^' 



(134) 



However, before application of this effective potential we must worry about the 
stricture of unitarity, which requires that 

2ImT = T^T (135) 

In the case of the S-wave partial wave amplitude mt{k) this condition reads 

lmmt{k) = k\mt{k)\'^ (136) 

and requires the form 



mt{k) = le'^'^'''> sin dt{k) (137) 



Since at zero energy we have 



lim mt{k) = —at (138) 
fe— »o 



It is clear that unitarity can be enforced by modifying this lowest order result 
via 

which is the lowest order effective range result. Equivalcntly, this can easily be 
derived in an effective field theory (EFT) formalism. In this case the lowest 
order contact interaction 

Tot = CM (140) 
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becomes, when summed to all orders in the scattering series, 

got(M) ^ M 1 ..... 



Identifying the scattering length via 

1 4tt 



at MCM 



(142) 



and noting the relation mt{k) = —^Tt{k) connecting the scattering and tran- 
sition matrices, we see that Eqs. 11391 and 11411 are identical. 

So far, so good. However, things become more interesting in the case of the 
parity-violating transitions. In this case the requirement of unitarity reads, e.g., 
for the case of scattering in the "^Si channel 

Imdt(fc) = k{m;{k)dtik) + d;{k)mp{k)) (143) 

where mp{k) is the ^Pi analog of the mt{k). Eq. 11431 is satisfied by the solution 

dt{k) = \dt{k)\e'^^'si^''^+^'^i^''^'> (144) 

i.e., the phase of the amplitude should be the sum of the strong interaction 
phases in the incoming and outgoing channels (6^. At very low energy we can 
neglect P-wave scattering, and can write 

ct{k) ~ ptmt{k), d;(fc) ~ A>i,(fc), dt{k) ~ Xtmtik) (145) 

This result is also easily seen in the language of EFT, wherein the full transi- 
tion matrix must include the weak amplitude to lowest order accompanied by 
rescattering in both incoming and outgoing channels to all orders in the strong 
interaction. If we represent the lowest order weak contact interaction as 

rotp(fc) = Dotp{n){ffi - a2) ■ {k + k') (146) 
then the full amplitude is given by 

^"'^ - (1 - c,Ma,.m'- cMom) - - 

where we have introduced a lowest order contact term C'op which describes the 
""^Pi-wave nn interaction. Since the phase of the combination 1 — Co(/i)G'o(fc) 
is simply the negative of the strong interaction phase the unitarity stricture is 
clear, and we can define the physical transition amplitude Atp via 



Dotpip) _ Atp 



(1 - Cot(/i)Go(fc))(l - Cop(M)Gi(fc)) (1 + ikat)il + ikHp) 



(148) 



Making the identification At = —^Atp and noting that 



^ cos5t(fc)e''^*('=) 



1 -I- ikat 
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then At is seen to be identical to the R-matrix element defined by Miller and 
DriscoUlSSl. 

Now that we have developed a fully unitary transition amplitude we can 
calculate observables. For simplicity we begin with nn scattering. In this case 
the Pauli principle demands that the initial state must be purely ^Sq at low 
energy. One can imagine longitudinally polarizing one of the neutrons and 
measuring the total scattering cross section on an unpolarized target. Since 
a • fc is odd under parity, the cross section can depend on the helicity only if 
parity is violated. Using trace techniques the helicity correlated cross section 
can easily be found. Since the initial state must be in a spin singlet we have 



|m,(fc)|2 ± AkRem*{k)d'^"{k) + 0{dl) (149) 



Defining the asymmetry via the sum and difference of such helicity cross sections 
and neglecting the tiny P-wave scattering, we have then 

A = ^ = 4fcAr (150) 

cr+ + (T_ 2|ms(fc)p 

Thus the helicity correlated nn-scattering asymmetry provides a direct measure 
of the parity- violating parameter A"" . Note that in the theoretical evaluation of 
the asymmetry, since the total cross section is involved some investigators opt 
to utilize the optical theorem via|f)4j.[()Kj 

Im ms(fc) 

which, using our unitarized forms, is completely equivalent to Eq. 11501 

Of course, nn-scattering is purely a gedanken experiment and we have dis- 
cussed it only as a warmup to the real problem — pp scattering, which introduces 
the complications associated with the Coulomb interaction. In spite of this com- 
plication, the calculation proceeds quite in parallel to the discussion above with 
obvious modifications. Specifically, as shown in |65) the unitarized the scattering 
amplitude now has the form 

Af Co.C,^(7y+(fc))exp2»cTo 
"^^^^^-'4^ l-C,,Gc{k) ^^^^^ 



where rj^{k) ~ Ala /2k and 



is the usual Sommerfeld factor and ctq = argr(£ -|- 1 -|- iri{k)) is the Coulomb 
phase shift. Of course, the free Green's function Go(fc) has also been replaced 
by its Coulomb analog 

f d^s C^iv+jk)) 
J (2^)'l7-i7+*£ 
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Remarkably this integral can be performed analytically and the result is 

Gcik) = [a^ + Ma (H{ir^+{k)) - log ^ - c)] (155) 
Here C is defined in terms of the Euler constant via C = 27r — 7^ and 

H(x)^il}{x) + — -\ogx (156) 
2x 

The resultant scattering amplitude has the form 

C2(,7+(fc))e2-o 



ms{k) 



ao 

where we have defined 

1 47r 



-life - M - Ma [H{r^+{k)) - log - C] 

C,^(r]+(fc))e^-° 

- Ma [h{ri+{k)) - log ^ - C] - ^kCl{ri+{k)) 



(157) 



aos MCqs 



fi, and h{r]+{k)) ^ReH{ir]+{k)) (158) 



The experimental scattering length aps in the presence of the Coulomb inter- 
action is defined via 

- — = -— + Ma (log ^-C) (159) 
acs aos V nA'la / 

in which case the scattering amplitude takes its traditional form 

C^(ri+(k))e^"'« 
--^,-MaH{n^j^{k)) 

Of course, this means that the Coulomb-corrected scattering length is different 
from its non-Coulomb analog, and comparison of the experimental pp scattering 
length — app = —7.82 fm — with its nn analog — a„„ = —18.8 fm — is roughly 
consistent with Eq. 11591 if a reasonable cutoff, say /j, ~ IGeV is chosen. Having 
unitarized the strong scattering amplitude, we can proceed similarly for its 
parity-violating analog. Again summing the rescattering bubbles and neglecting 
the small p-wave scattering, we find for the unitarized weak amplitude 

_ iJo,p(M)fi^(77+(fc))e'(-°+-^) _ Ac,pC^^(y?+(fc))e»(-°+-^) 

\-CMGc{k) ~ -^-MaH{irj+ik)) ^ ^ 

Here again, the DriscoU-Miller procedure identifies Acsp = Af^ via the R-matrix. 
Having obtained fully unitarized forms, we can then proceed to evaluate the 
helicity correlated cross sections, finding as before 

A. ^ ^-±^ ^ 8fcRe<(fcKnfe) ^ ,,^p (,g2) 
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Note here that the superscript pp has been added, in order account for the 
feature that in the presence of Coulomb interactions the parity mixing parameter 
As which is appropriate for neutral scattering is modified, in much the same way 
as the scattering length in the pp channel is modified (c/. Eq. I159|) . On the 
experimental side such asymmetries have been measured both at low energy 
(13.6 and 45 MeV) as well as at higher energy (221 and 800 MeV) but it is only 
the low energy results^ 

A,,(13.6MeV) = -(0.93 ± 0.20 ± 0.05) x lO^^HT] 

(45 MeV) = -(1.57 ± 0.23) X lO-^'EHI (164) 

which are appropriate for our analysis. Note that one consistency check on these 
results is that if the simple discussion given above is correct the two numbers 
should be approximately related by the kinematic factor^ 

A(45 MeV)/A(13.6 MeV) ~ /ci/fcz = 1.8 (165) 

and the quoted numbers are quite consistent with this requirement. We can 
then extract the experimental number for the singlet mixing parameter as 

= ^ = -(4.0 ± 0.8) X 10-^ fm (166) 
4fc 

In principle one could extract the triplet parameters by a careful nd scattering 
measurement. However, extraction of the neeeded np amplitude involves a de- 
tailed theoretical analysis which has yet not been performed. Thus instead we 
discuss the case of electromagnetic interactions and consider np ^ dj. 



7 Parity Violating Electromagnetic Interaction: 

np ^ dj 

A second important low energy probe of hadronic parity violation can be found 
within the electromagnetic transtion np ^ dj. Here the np scattering states 
include both spin-singlet and -triplet components and we must include a bound 
state — the deuteron. Analysis of the corresponding parity-conserving situation 
has been given previously, so we concentrate here on the parity violating situ- 
ation. In this case, the mixing of the scattering states has already been given 
in Eq. 11301 while for the deuteron the result can be found from demanding 
orthogonality with the "^^i scattering state — 



Mr) = fl + Pt{<^p + >?«) • -iV -t- Xtiap - a„) • J -^-e-^"- (167) 



*Note that the 13.6 MeV Bonn measurement is fully consistent with the earlier but less 
precise number 

Aft = -(1.7 ± 0.8) X lO-'^EHl (163) 

determined at LANL. 

^There is an additional fc-dependence arising from As but this is small. 
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Having found Ag via the pp scattering asymmetry, we now need to focus on the 
determination of the parity- violating triplet parameters pt , AJ . In order to do 
so, we must evaluate new matrix elements. There are in general two types of 
PV El matrix elements, which we can write as 



Amp = (HeixI{^p - ^n)Xt + SeixK^p + ^n)Xt) • Co 



(168) 



and there exist two separate contributions to each of these amplitudes, depend- 
ing upon whether the parity mixing occurs in the inital or final state. We begin 
with the matrix element which connects the ^Pi admixture of the deuteron to 
the ^5o scattering state. 



HeiCP-'S) = 



espXt 1 
2 3 
eso47rAt 



^ io 



(r)f • Vipd{r) 



6 



esQ^/SiT'yXt e'' 



drr^ — (sin kr cos Sg + cos kr sin 5^ )r • V 
kr 



k 



X (1 — 7-;— ) / dr(sin fcrcos^s + cosfcrsin^s)e 

dl Jo 

esoAtVSTFye"'*^""''^"'' 



6^1 + k^al 



(fc2+72) 



2 - Ids 



(fc2+72)2 



Equivalently we can use EFT methods using the diagrams of Figure 4. We have 
from Figure 4a 



espXt 1 
2 3 



d^q 1(0)* 



(27r)3 ^fe 

6 (fc2+72)2 

while from the bubble sum in Figure 4b-|-c we find 



(170) 



espXt 1 
2 3 
esoM^ySn^Xt 



C'os 



d\ 

l-CosG„{k)JJ (27r)3 



Q 



Os 



d^q 



J (27r 



1 - CosGoik) 
9^ + 372 



0' (q'+l'n'i-^+ie) 
Here the integral may be evaluated via 



Go(f=0,g)V,~-(#4g)) 



(171) 



/ 



d^q 



g2 + 37' 



(27r)3 (g2+^2)2(fei_£+i,) 



= (1-27^37^) 



1 



d^2 ^ (/,2 + ^2) 



(Go(fc)-Go(z7)) 
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47r(7-iA;)2 ^ ' 

Summing the two results we find 



2 3(fc2+72)2\^ ' 1 + ikas 

6 (fc2 + 72)2(1 +ifca,) 



(173) 



as found using coordinate space methods. 

The matrix element Hei also receives contributions from the El amplitude 
connecting the deuteron wavefunction with the mixture of the final state 
wavefunction mixed into the ^Sq. This admixture can be read off from the 
Green's function representation of the scattering amplitude as 

Sspipiso = -ims{k){ap - an) ■ (l'^4) 
and leads to an El amplitude 

HEr{'P-'S) = ££2^^|^i|d3,V.(r)f.V^^r,„(fcr) 



egoA7VB7T7 

6 V 1 ~ ^^^s 



/ dr{l + ikr) 
Jo 



g (7+ifc)r 



6 \1 — ikas J (7 + ifc)2 



itan ^kas 

(7(72 + 3fc2) - 2iA;3) (175) 



6(fc2+ 72)2^1 + /c2a2 

Equivalently we can use EFT techniques. In this case there is no analog of 
Figure 4a. For the remaining diagrams, however, we find 



2 i \l - CosGo{k) 
(2^) 



/ (0^c;(?)V,-gGo(r = O,g-) 



1 - Co«Go(fc) 

- ( 

esoA7 / Cos \ M 7(72 + 3A;2) - 2ik'^ 



l-CosGo{k)J A-K i^^ + k^y 



(176) 



I.e., 
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as found in coordinate space. 

The full matrix element is then found by combining the singlet and triplet 
mixing contributions — 

Hei = HEieP-''S) + HEiCP-'S) 

6^1 + Pa2(fc2+ 7^)2 
X [At(fc2 + 372-2a,73) + A^P7a,[(72 + 3fc2)7-2iA:3])] (178) 

In the case of the PV El matrix element Sei the calculation appears to be 
nearly identical, except for the feature that now the spin triplet final state is 
involved, so that the calculation already performed in the case of Hei can be 
taken over directly provided that we make the substitutions A"^, Xt pt, as — > 
at- The result is found then to be 

Se, = "7 T^2n2 ^ [(^' + - + «*[7(7^ + 3fc^) - 2*^^])] 

(179) 

At the level of approximation we are working we can identify at with I/7 so 
that Eq. 11791 becomes 

-1 k 



e.o^/8^e-*-" -p,(72 + 2fc2-^^) 

JEi = , (18U) 

3^1 + ^(fc2+ 7^)2 

Now consider how to detect these PV amplitudes. The parity violating 
electric dipole amplitude Hei can be measured by looking at the circular po- 
larization which results from unpolarized radiative capture at threshold or by 
the asymmetry resulting from the scattering of polarized photons in photodis- 
integration. At threshold, we have for photons of positive/negative helicity 

Amp_|_ = {±Gmiv + HEi)e^ ■ xl{^p - $n)xt + Ssii^ ■ xK^p + ^n)xt (181) 
and the corresponding cross sections are found to be 

a±{np^dj) = --^^\^Gmiv+Hei\^ + 0{SI,) (182) 

Thus the spin-conserving El amplitude Sei does not interfere with the leading 
Ml and the circular polarization is given by 

P-i = T = = ~^ T\ \(h2^ 2\ V^iV^ + 37 - 2as7 ) 

+ A,"f7a.(7'+3fc'))] (183) 

A bit of thought makes it clear that this is also the asymmetry parameter 
between right- and left-handed circularly polarized cross sections in the photo- 
disintegration reaction 7±(i — s- np, and so we have the usual identity between 
polarization and asymmetry which is guaranteed by time reversal invariance 



P-y {np d'y) ~ Ay (-fd np) 
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In order to gain sensitivity to the matrix element Sei one must use polarized 
neutrons. In this case the appropriate trace is found to be 



Tri(3 + ai ■ a2)T^{l + ■ n)T^ = A\Gmiv\'' + 8Re Qnv^Bin • (184) 

In this case Hei does not interfere with Gmiv and the corresponding front-back 
photon asymmetry is 

^ 2ReGl,^ySEi ^ 8Mpt + 2^ 



In principle then precise experiments measuring the circular polarization and 
photon asymmetry in thermal neutron capture on protons can produce the 
remaining two low energy parity violating parameters A"'' and pt which we 
seek. At the present time only upper limits exist, however. In the case of the 
circular polarization we have the number from a Gatchina measurement [75] 

= (1.8± 1.8) X 10"^ (186) 

while in the case of the asymmetry we have 

= (-1.5±4.7) X 10"^ (187) 

from a Grenoble experiment j71|. This situation should soon change, as a new 
high precision asymmetry measurement at LANL is being run which seeks to 
improve the previous precision by an order of magnitude|72|. 

Appendix 

Of course, as we move above threshold we must also include the parity- 
violating Ml matrix elements, which interfere with the leading El amplitude 
and are of two types. The first is the Ml amplitude which connects the ^Pi 
admixture of the deuteron with the np scattering state as well as the Ml 
amplitude connecting the ^5*0 admixture of the final '^P scattering state with 
the deuteron ground state. For the former we have^ 



Amp = Jlnxli^i - ^2) ■ £7 X s^{ai - ^2) • kxt (188) 



where 



while for the latter we find 

Amp = JmixK^i - ^2) ■ k{ai - CT2) • £7 x s^xt (190) 

^For simplicity here we include only the dominant isovector Ml amplitude. A complete 
discussion should also include the corresponding isoscalar Ml transition. 
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where 



t'' = 

Ml 4JV^ "-avvy - ■ rav/ (1 - ifca,) (7 + 



(191) 

A second category of PV Ml amplitudes involves that which connects the 
^Pi piece of the deuteron wavefunction with the ^Pi or ^P np scattering states 
as well as the Ml amplitude connecting the or admixture of the final 
^P scattering state with the deuteron ground state. For the former we have 

KmixIH^i - ^2) ■ £7 X ^7(^1 + ^2) • kXt (192) 

^MixlWi + ^2) • £7 X ^7(^1 + ^2) ■ kxt (193) 



where 



while for the latter we have 



Amp = Kl^ixli^i + CT2) • k{^i - ^2) ■ £7 X s^xt (196) 



Amp = iMiXt (cti + 0^2) • A;(cti + CT2) • Cj x s^xt (197) 



where 



= y 'r' "^'^'^ = - AM (l-^fca,)(7 + ^fc) ^'''^ 

iMi = ^P*m,(fc)y dr-e W = - (1 _ + (199) 

To this order we can use at — I/7, so that 

^^1- 4M fc2 + ^2 (200) 

jb _ esQps PtVS^ . . 

4M fc2+^2 (201) 

We see then that this piece of the Ml amplitude has the form 

Amp = 2i(e% x s^) x • X/[^mi(^i - ^2) + iMi(^i + '?2)]Xt (202) 

The relevant traces here are 

Tri(3 + (71 • CT2) (Jmi('?i - ^2) ■ £7 X Sj{ai - 0^2) • 
+ Jliii^i - ^2) ■ k{ai - (72) • £7 X s^^ i(l + 0^2 • n) 
~ 2(>/jvfi + '^Mi)£7 X • A; + lii^Jj^-^ Jmi)'^ ' (£7 x S7) x k (203) 
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Tr^(3 + (Ji • (72) (KMii^i - CT2) • £7 X s^(cti + S2) -k 

+ K\,i^{ai + (72) • ^(cti - CT2) • £7 X s^^ -(1 + (72 • n) 

= 2(L^]^ + L\ji)t^ X Sy ■ k ~\- 2z(L^,j]^ — L\^i)n ■ (e-y x s^) x fc 

= AiLlj^n ■ (e^ x s^) x fc (204) 

Tri(3 + (7i • $2) (L%ji{ai + G2) ■ x s~^{gx + (72) • fc 



4^ 



1 



+ -LAfi((7i + 02) ■ k{ai + (72) • £7 X s-yj -(1 + (72 • n) 

— '^(^■^Mi ^" ^ Mi)^i X s-y • /c -t" 'lii^L ^j-^ L^j-y^fi • (c-y X s^) X k 
= 4:iL%f^n ■ {ij x s^) x k (205) 

and the corresponding contribution to the cross section for photo(iisintegration 
by photons of differing helicity is 

"""^ " 1271^72 + fc2) (l*^^!'^ =^ ^^(^G^iWii + An) 

- 3(fc2+72)3 ^ 9(^ + 7^) I * ~ ^ ""^TTP^ j ^^^^^ 

Note that there is only sensitivity to the couphngs Jmi here. In order to 
have sensitivity to the couphngs Kmi, Lmi we must look at the El contribution 
to the cross section for the radiative capture of polarized neutrons- 

^ = iS\vrl\M - - + 8ReG^i(j^jtfi + ^Mi)g7 ■ ks-, xk-n) 

3a7 k--{k- s,f , 2«7(Mv + ^^s)pt_.^ . x fc • n (207) 



M3|zvei| fc2 + 72 M^\n 



rel I 



Careful analysis of above threshold experiments should include these NLO 
corrections to the analysis. 

8 Connecting with Theory 

Having a form of the weak parity- violating potential V^^^{r) it is, of course, 
essential to complete the process by connecting with the S- matrix — i.e., ex- 
pressing the phenomenological parameters , pt defined in Eq. 11341 in terms of 
the fundamental ones — C^, Ci defined in Eq. ^ This is a major undertaking 
and should involve the latest and best NN wavefunctions such as Argonne V18. 
The work is underway, but it will be some time until this process is completed. 
Even after this connection has been completed, the results will be numerical 
in form. However, it is very useful to have an analytic form by which to un- 
derstand the basic physics of this transformation and by which to make simple 
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numerical estimates. For this purpose we shall employ simple phenomenological 
NN wavefunctions, as described below. 

Examination of the scattering matrix Eq. 11301 reveals that the parameters 
As,t are associated with the (short-distance) component while pt contains con- 
tributions from the both (long-distance) pion exchange as well as short distance 
effects. In the former case, since the interaction is short ranged we can use this 
feature in order to simplify the analysis. Thus, we can determine the shift in the 
deuteron wavefunction associated with parity violation by demanding orthogo- 
nality with the ^Si scattering state, which yields, using the simple asymptotic 
form of the bound state wavefunction|42 | . |43 | 

Mr)= [l + /9t(CTp + CT„)--iV + At(ap-a„)--zV)l J —e-''- (208) 

J V ZTT r 

where = 2.23 MeV is the deuteron binding energy. Now the shift gener- 

ated by y^^(r) is found to be|ll],|12I 



SMr) - I d\'G[r,r^)V''''{7^)Mr') 



(^-^^ ■ J d3/r'F^^(r')^4r') (209) 



where the last step is permitted by the short range of V^'^ (f). Comparing Eqs. 
12091 and 12081 vields then the identification 

^ X,Xt ^ I d'r'ia, - a^) ■ r'V'''' {r^)Mr')Xt^o (210) 



2tt Wit 

where we have included the normalized isospin wavefunction ^0 since the po- 
tential involves ti , T2 ■ When operating on such an isosinglet np state the PV 
potential can be written as 

V'^'ir') = [(Ci-3C3)(ai-a2)-(-»V/™(r)+2/„,(r).-iV) 

X 

+ (Ci - 3C'3)(ai X as) • V/„,(r)] (211) 
where fm{r) is the Yukawa form 

Mr) 



defined in Eq. Using the identity 



((?! X (T2)i(l +<Ti • (Ta) = i(CTl - CT2) (212) 
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Eq. l21Ul becomes 



7 ^ 2M An 



167rA3 3 



-2(3C3 - Ci)/„(r)^%^^ + (3(73 - 5Cs - C, + C^)^^Mr) 



dr 



dr 



7 1 2Mm2 / 2to(6C3 - SC^ - 2Ci + Ci) + 7(15C3 - 30^ - 5Ci + Ci) 

2^' ^*12 47rA3 1^ (7 + m)^ 



or 



(213) 



A* 



Mm2 / 2m(6C3 - 3(73 - 2Ci + (7i) + 7(15(73 - 3(73 - 5(7i + (7i) 
67rA| 



(7 + m)2 



(214) 

In order to determine the singlet parameter A"^, we must use the ^Sq np- 
scattering wavefunction instead of the deuteron, but the procedure is similar, 
yieldinggg,!!! 



d7{k)xs ^ I d\'{a, - CJ2) ■ r'^^^(r')^.s,(r')x.a 

and we can proceed similarly. In this case the potential becomes 



(215) 



A3 

X 



[Ci + C3 + \c^){Si - (72) • (-^V/™(r) + 2/„(r) • -iV) 



1 ^ 



+ (Ci+C73 + 7TC5)(ai xa2)-V/™(r) 
6 



(216) 



and Eq. 12 151 is found to have the form 



X 2(Ci+C3+4C5)/„.(r) 



dr 



+ (Ci+Ci+C3 + C3+4(C5 + C5))^^^i5oW 



1 2Mto2 



36 47rA3 



1 



(fc2 + to2)2 

X (^cos(5,(4fc2(Ci + C3 + 4(75) + [Ci +Ci+C^ + Cz+ 4((75 + C^)){k'' + 3™^)) 

2r77 _ _ 

+ — sin5,(((7i + C3 + 4C5)(m2 + ik^) + [d + (7i + C3 + (73 + 4((75 + Cs))™^) 



(217) 
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which, in the hmit as fc ^ 0, yields the predicted value for A"^- 
1 .. M 



as k~*0 DTTfls A?, 



3(Ci +C1+C3 + C3+ 4(C5 + C5)) 



2ma^P(2C7i + Ci + 2C3 + C3 + 4(2^5 + C5))) 



(218) 



Similarly, we may identify 
1 



M 



3(Ci + Ci + C2 + C2 + C3 + C3 + C4 + C4 - 2(C5 + C5)) 
2maPP(2Ci + Ci + 2C72 + C2 + 2C3 + C3 + 2C4 + C4 - 2(2^5 + C5))) 

3(Ci + Ci - C2 - C2 + C3 + - C4 - 6-4 - 2(^5 + C5)) 



1 M 
lim dr(fc) = 7 TT 



2ma™(2Ci + Ci - 2C72 - C2 + 2C3 + C3 - 2C4 - C4 ^ 2{2C5 + C5))) 



In order to evaluate the spin-conserving amplitude pt, we shall assume domi- 
nance of the long range pion component. The shift in the deuteron wavefunction 
is given by 



but now with'' 



a / dV'Go(f,r-)T//^(r')^,(r') 



M 
47r " 



h-^gTrNN 1 

V2Mml 2 



~7|r— r 



(220) 



(222) 



Of course, the meson which is exchanged is the pion so the short range as- 
sumption which permitted the replacement in Eq. I2U9I is not valid and we must 
perform the integration exactly. This process is straightforward but tedious ^7]. 
Nevertheless, we can get a rough estimate by making a "heavy pion" approxi- 
mation, whereby we can identify the constant pt via 



27r 

which leads to [HI 
7 

PtXt 



PtXt 



. M 
'32^ 



(223) 



^Here we have used the identity 



1 47r , ^ n2 ^Kg-nNN 



dr 



(219) 



(n X T2) = - + ti ■ fa) 



(221) 
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fjlfi 1 Kg^NN 7 + 2m^ 
~ V 277^^*9671 72 {l + m^Y ^^^^^ 

We find then the prediction 

pi = ^7 — ; rj/iir (225) 

At this point it is useful to obtain rough numerical estimates. This can 
be done by use of the numerical estimates given in Table 2. To make things 
tractable, we shall use the best values given therein. Since we are after only 
rough estimates and since the best values assume the DDH relationship — Eq. 
IHl between the tilde- and non-tilde- quantities, we shall express our results in 
terms of only the non-tilde numbers. Of course, a future complete evaluation 
should include the full dependence. Of course, these predictions are only within 
a model, but they has the advantage of allowing connection with previous the- 
oretical estimates. In this way, we find the predictions 

\t = [-O.O92C3 - O.OWCi] 

= [-0.087(C3 + 4C5) - 0.037Ci] m'^ 

Af = [-0.087(C3 + C4-2C5)- 0.037(Ci-)-C2)]m-i 

XT ^ [-0.087(C3 - C4 - 2C5) - 0.037(Ci - C2)] m'^ 

pt = Q.MQKm:^'^ (226) 

so that, using the relations Eq. Inland the best values from Table 1 we estimate 



(227) 
(228) 

At this point we note, however, that is an order of magnitude larger than 
the experimentally determined number, Eq. 11661 The problem here is not with 
the couplings but with an important piece of physics which has thus far been 
neglected — short distance effects. There are two issues here. One is that the 
deuteron and NN wavefunctions should be modified at short distances from the 
simple asymptotic form used up until this point in order to account for finite 
size effects. The second is the well-known feature of the Jastrow correlations 
that suppress the nucleon-nucleon wavefunction at short distance. 

In order to deal approximately with the short distance properties of the 
deuteron wavefunction, we modify the exponential form to become constant 
inside the deuteron radius i?.j42 | .pH ] 

X_e-7-^iV "^-^ (229) 

2tt r \ -e ' r > R 



At - 


-2.39 X 10" 




-3.41 X 10"'^ fm 


= 


-1.12 X 10" 




-1.60 X 10"^ fm 


= 


-3.58 X 10" 




-5.22 X lO"'^ fm 


A?" = 


-2.97 X 10" 




-4.33 X lO"'^ fm 


Pt = 


1.50 X 10"^ 


= 2. 


.14 X lO"'^ fm 



42 



where 

^ I 7 exp7i? 



is the modified normaUzation factor and we use R=1.6 fm. For the NN wave- 
function wc use 



= { . ^P^'- " ' (230) 



pr -^+ip 



r > r. 



where we choose = 2.73 fm and pors = 1.5. The normaUzation constant A{p) 
is found by requiring continuity of the wavefunction and its first derivative at 

A(p) = Vp'^ +Po'^'> sin pr s - pg^ cos pr s 
sin p2^^ prsil + ipas ) 

As to the Jastrow correlations we multiply the wavefunction by the simple phe- 
nomenological formj45j 

0(r) = 1 - ce"'*''', with c = 0.6, d = 3 fm"^ (232) 

With these modifications we find the much more reasonable values for the con- 
stants XPP "P and At 

XPP = [-0.011(C3 + C'4-2C5)- 0.004(Ci-fC2)]m-i 

XT = [-0.011(^3 -C4 + 2C5)- 0.004(Ci-C2)]m;i 

= [-0.011(C3 + 4C5) - 0.004Ci] 

At = [-O.OI9C3 - 0.0003Ci]to-i (233) 

Using the best values from Table 2 we find then the benchmark values 



XPP - 


-4.2 X 10" 




-6.1 X 10" 


-8 


fm 


Ar = 


-3.6 X 10" 




-5.3 X 10" 


-8 


fm 


A«f - 


-1.3 X 10" 


-v-i = 


-1.9 X 10" 


-8 


fm 


At = 


-4.7 X 10" 




-6.7 X 10" 


-8 


fm 



(234) 

Since pt is a long distance effect, we use the same value as calculated previously 
as our benchmark number 

Pt = 1.50 X lO"'^ = 2.14 x 10""^ fm (235) 

Obviously the value of Af^ is now in much better agreement with the ex- 
perimental value Eq. 11661 Of course, our rough estimate is no substitute for a 
reliable state of the art wavefunction evaluation. This has been done recently 
by Carlson et al. and yields, using the Argonne V18 wave functions |7Hj 

XPP = [-0.008(C3 + C4 - 2C5) - 0.003(Ci + C2)] (236) 
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in reasonable agreement with the value calculated in Eq. 12331 Similar efforts 
should be directed toward evaluation of the remaining parameters using the best 
modern wavefunctions. 

We end our brief discussion here, but clearly this was merely a simplistic 
model calculation. It is important to complete this process by using the best 
contemporary nucleon-nucleon wavefunctions with the most general EFT po- 
tential developed above, in order to allow the best possible restrictions to be 
placed on the unknown counterterms. 

9 Summary 

For nearly fifty years both theorists and experimentalists have been struggling 
to obtain an understanding of the parity-violating nucleon-nucleon interaction 
and its manifestations in hadronic parity violation. Despite a great deal of effort 
on both fronts, at the present time there still exists a great deal of confusion 
both as to whether the DDH picture is able to explain the data which exists and 
even if this is the case as to the size of the basic weak couplings. For this reason 
it has recently been advocated to employ and effective field theory approach 
to the low energy data, which must be describable in terms of five elementary 
S-=matrix elements. Above we have discussed both connection between these S- 
matrix elements and observables in the NN system via simple analytic methods 
based both on a conventional wavefunction approach as well as on effective field 
theory methods. While the results are only approximate, they are in reasonable 
agreement with those obtained via precision state of the art nonrelativistic po- 
tential calculations and serve we hope to aid in the understanding of the basic 
physics of the parity-violating NN sector. In this way it is hoped that the round 
of experiments which is currently underway can be used to produce a reliable 
set of weak couplings, which can in turn be used both in order to connect with 
more fundamental theory such as QCD as well as to provide a solid basis for 
calculations wherein such hadronic parity violation is acting. 
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